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Preface 


This document records a simple project which I and one of my classmates took on over 
several weeks of our Winter Break, just on a whim. We’re both Math/Physics concentra- 
tors, but we take a special interest in potentially transformative technologies. Elon Musk’s 
Hyper loop transportation system, with its promise for cheap, fast, and green transportation, 
falls into that category. For more details on the proposal, check it out online via Google 
Search. 

There were many aspects of the Hyperloop project that we could have worked on, but 
I felt that route optimization was an area that could play to our prior experience with 
functional analysis on abstract spaces, a course which both of us thoroughly enjoyed, be- 
cause, I must say a cliche, the astonishingly powerful and beatiful ideas expanded our minds. 

Therefore our journey begins at the most abstract aspect of the Hyperloop project: de- 
termining the route. What did Musk have to say about this? Well, in his Alpha document 
describing the Hyperloop proposal, Elon Musk stated that he was looking for a trip path 
and velocity that did not exceed .5 g’s at any point. Our goal was to minimize trip time 
and cost under this passenger comfort constraint, using the path and velocity profile as our 
variables. 

The results are, to say the least, interesting. Over the course of just a few weeks, Jon 
and I learned an incredible amount in such diverse topics, from computer science, to math- 
ematics, to physics. We both hope you find beauty as well, as you navigate this document. 

Enjoy! 
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Chapter 1 

Calculating Time T('y) 


1.1 The G - force Constraint and Resulting Equation 


Presume a proposed Hyperloop path 7 : M — y M 3 is parametrized by the distance along the 
path. We can then express 7 in an Earth-Centered, Earth-Fixed coordinate system: 


7(7) = ( x(s),y(s),z(s )) 


Then we have that the velocity profile v(s), bounded by the maximum centripetal accelera- 
tion tolerance g, must satisfy the following inequality, with g rad and .g ]in being the centripetal 
and linear accelerations respectively: 


g<\ 9ia,d + f hi 


( dv(s ) v( sj“\ 

(fl'lin ~dt~ ’ 9rad ~ ) 


_ v(s) 2 


To hnd the optimal velocity prohle we take the effective acceleration to be the maximum 
allowed value, which is the constant G, so that we have the equality case in equation (1). 

2 

G 2 = 9 rad + 9hn = 


(s = vt) 



Therefore our g-force constraint yields a differential equation, whose solution, v(s), is the 
optimal velocity prohle we are looking for. The radius of curvature is a function of the path 
in the following way: 


r(s) 


(x' 2 + y' 2 + z' 2 ) I 

a/ {z"y' — y"z') 2 + (x"z' — z"x') 2 + (y"x' — x"y') 2 


1.2 Solving the Equation Locally 

The boxed differential equation of the previous section is incredibly difficult to solve, because 
r(s) is completely arbitrary, even chaotic. We therefore solve the simplest case, where r(s) 
is constant, and break up the path into a series of regions of constant curvature. 
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1.2.1 Calculating Radii of Curvature 

The first step, of course, is to convert a path into a series of sections of constant curvature. 
The process has three steps, and involves the internet, Google Earth, and Microsoft Excel: 

(a) In Google Earth, get driving directions from point A to point B (This will be our path 
7 ), and export the route as a .kml hie. 

(b) Visit http : //www. geocontext . org/publ/2010/04/prof iler/en/?import=kmV and 
import your .kml hie, and click the .csv option: 


I 


I 




Q Help 

How to make a topographic profile? 


1. Reset 

2. Find your area of interest on the map 

3. Select the cursor min. 2 points (max. 300) 

4. Ready - site profile will be generated in 
seconds 

a. Embed the chart on your site 

b. Copy and save the link to the chart 

c. Add the route to the map 

Program Geocontext-Profiler allows you to make 
topographic profiles anywhere on Earth in the 
seabed and ocean floor. Can be widely used: in 
the natural sciences (geomorphology, 
hydrography) for the education and tourism - 
hiking route planning, bicycle, car. 

Within the program, you can find some advanced 
options that allow you to create a profile along the 
road, bicycle and pedestrian paths, and measuring 
the slope angle. The program can import KML, 
KMZ and GPXfrom GPS devices. 

For educational purposes several pre-programmed 
profiles of interesting geographic features, such 


(c) Copy and paste the results into Microsoft Excel, and after pasting, click the clipboard 
for paste options, to allow Excel to recognize the .csv format of the contents. Now save 
the file as a .csv. 

(d) This gives us a series of longitude and lattitude points with corresponding elevation 
(meters above sea level) {(9 i: fa, ej)}” =1 . These coordinates denote points on our path 
7 . We then transform into an Earth-Centered, Earth-Fixed coordinate system, in order 
to calculate curvature. 

Yi = ( Xi,yi,Zi ) = (R® + ei) (cos 9 i cos cos 9 i sin fa, sin 9 i) G M 3 

(i?© = Radius of the Earth) 

(e) We then approximate the radius of curvature R(ri) at the point r* by the radius of the 
circle containing rj_ 1; r*, r i+1 . This is a nonlinear system of three distance equations 
in 3D-space, 


R(ri) = ||r;_i - c 1 1 = ||r< - c|| = ||r m - c|| 

where c = (x, y, z) is the position of the center of the circle. There is now one last step: 
we want to describe the path as a series of regions of constant curvature, but we have 
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defined curvature for only the points. So for the interval [sj,s, + i], we assign a radius 
of curvature as follow^] 


7?( [sj, ) 


i — 1 or n — 1 

i?(rj_i) + R(vi)\/2 otherwise 


1.2.2 Solving the Equation with Constant Curvature 


We have made the discrete approximation of r(s) precisely so that we can solve our differential 
equation with r(s) = R, a constant. Our solution then describes the optimal speed of the 
pod entering a circular segment with intitial velocity Vq'- 


g 2 = + „ (s y- 

1 ( du 

“il * 1 + 


dv(s] 


ds 


f du 
\ds 


= AG 2 


ds= 2G 


R 
~ ~2 
R 

s = — Sill 
2 


u(s)' 2 

u(s) 

RG 

du 


+ C 


1 - 

dx 


\/l — 


R 2 G 2 

= + C 


x- 


-i 


u 


+ C 


uis ) = RG sin 


RG 


(u = v 2 ) 

(Roberts- Ward Equation) 


(x = u/ RG) 


We now solve for C based on the initial condition w(0) = uq\ 


(2 C 


u 0 = —RG sin ( — ) — » C = — sin 1 


V R 

u(s) — RG sin 


R 


up \ 
RGJ 


2 s . _j / Up \ 

h Sill 

R \RGJ 



v(s) = y / RG sin 

2 s . _i f v% \ 

R +Sm \RG J 



(Roberts- Ward Solution) 


Note that this equation does not account for a maximum velocity prescribed by any other 
considerations such as thermodynamic and stress considerations, shockwaves. Furthermore 
this equation assumes that the velocity can be modified at will, in particular that there are 
no restrictions on the placement of the linear accelerators and also that the acceleration is 
smooth. 

^^Note that this formula is symmetric, forwards and backwards. This will turn out to be immensely useful 
when solving the differential equation, because the physics of the solution contains this symmetry. In other 
words, the path velocity should be the same whether you go from start to finish (forwards), or from finish 
to start (backwards). 


1.3 Solving the Equation Globally 

The goal now is to stitch together local solutions (for v(s)) to form a global optimal velocity 
curve. The idea is to break up this job into two steps. We first solve for the maximum values 
of the velocity at the boundaries of each segment of the path, i.e. the entering and exiting 
velocities for each section. Note that we have to be careful, for if such velocities are too far 
apart, and the section is too short, the resulting solution will be discontinuous. 


1.3.1 Finding Optimal Boundary Conditions 

(a) Focus on u; in particular, focus on {u(sj)}f =1 C M, i.e. the values of u at the boundaries. 
Once these boundary conditions are solved for (which is the goal of this algorithm), 
the analysis of the next subsection should be able to determine u(s) for all sGl. 

(b) Set {«(sj)}” =1 initially to be {G.R(sj)}" =1 . Of course there does not exist a u satisfying 
these n boundary conditions, but we have to have a starting point that aims high. In 
particular, we know that the actual solution to u has boundary values less than these 
initial ones (unless the entire route has uniform curvature). 


(c) We know that the algorithm has finished, (i.e. that a continuous u satisfying the 
boundary conditions exists) once, if for all i from 1 to n: 


u i+ 1 e [ui + A + {m ) , Ui - A_ (u^ 


(Notation: Ui = u(si )) 


Where A + (wj), A_(wj) are respectively the maximum positive and negative changes in 
u permitted clue to the physics of the region as a function of the initial square velocity 
Ui . More specifically, 


A + {v,i) = sup ( RiG sin 

.SG[0,As;] 


2s 

~R ~ Sm 


-l 


Ui 


RiG 


— u. 


A _(ui) = 


inf ( i?,;Gsin 

se[0,Asj] 


2s 

~R 


sm 


-l 


Ui 


RiG 


u, 


Where if A _(ttj) > iq, as an extra step we set A -(ui) = Ui, so as to not permit negative 
velocities and thus wasted time. 


(d) Our algorithm thus is the following: if the algorithm is not finished, then there exists 
some i such that 

ut+i £ [«i - A _(«»), + A + («i)] 

Now pick the least such i. Then execute the following: 


If 

Do 

Ui+ 1 > Ui + A + (ui) 

u i+ 1 < Ui - A _(ui) 

Set u i+ 1 = Ui + A+(wj) 
Set Ui = u i+ 1 + A _(wj) 


To see this algorithm in action, let’s now focus on what happens to the first 100 regions, 
for detail. Notice how the algorithm is making incredibly minor changes by the 3rd run 
through: 
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For a driving route from SF to LA, the boundary conditions converge after three iterations 
of the algorithm. 


1.3.2 Revisiting the Case of Constant Curvature 

Consider a region of curvature Rj lasting from s* to s 8+ i . By the above algorithm, we deter- 
mine the boundary conditions to be: u(si) = Ui, and w(sj + i) = tq + ±. Our goal in this section 
is to find the optimal curve u(s) connecting these two points. 


According to the Roberts- Ward solution the transition u- value between a boundary point 
u\ and a boundary point u -2 where s zero is the s value of the transition point between u\ and 
« 2 , and G is a parametric value for the maximum allowed acceleration is given by 


SolFun(ui, U 2 , s zero , G, s ) = 112 sin 


2 G u 2 

— (\s - s zero ) + — arcsm 

U 2 ZLr 



(1.1) 


A physically realistic solution requires that the pod steadily maintain its maximum allowed 
u- value after achieving it. The maximum allowed u- value is attained when the argument of 
the sin function in u(s) equals | 


2 G 

U 2 


s , U 2 . ( Ul 

•‘’max "Szero) T .. ^ t arcsill I 

ZLt \U 2 

x , u 2 . fu 1 

.S'niax "Szero) T .. ^ < arcsill I 

ZLt \U2 


3 max °zero/ 


3 max °zero/ 


"Smax 


7 T 
2 


U 2 7T 

2G2 


U 2 7T 

2G2 


u 2 

2 G 


arcsm 



U2_ 
2 G 



U2_ 

2 G 


( n 

arcsm 

V 2 



+ s 


zero 
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Now we can define a piecewise transition function, TransFun. 


There a three cases for TransFun to handle: U\ > u 2 , u 2 > Ui, and u 2 == U\. In the 
trivial case U\ == u 2 , the solution function is identically zero and Transfun returns U\ as 
the u- value. 

In the case where u \ < u 2 we can make Transfun(«i, u 2 , s zero , G, s) physical over the en- 
tire domain by setting the u-value equal to u\ in the region with .s < s zero , equal to 
SolFun(M!, u 2 , s zero , G, s) in the region with s zero < s < s max , and equal to u 2 in the re- 
gion with s max < s. 


In the case that that u 2 < u \ we can use the Roberts- Ward solution again by exploiting 
the symmetry of the physical system. If we reflect the boundary conditions across the line 
s — 0, then we obtain boundary conditions identical those in the case that u\ < u 2 . So by 
swapping the roles of u\ and u 2 in our solution function and reflecting s zero and s across the 
s = 0 line. So we make Transfun(rti, u 2 , s zero , G, s) physical by setting the u-value equal to 
u i in the region with s < —s m ax, equal to SolFun(rt 2 , Ui, — s zero , G, *^s) in the region with 
Szero > s > —s max , and equal to u 2 in the region with s zero < s. 


The next step is to define ConnectFun, which defines the u-value over a region with 
radius of curvature R and with boundary conditions u \ and u 2 . Given the myriad of detailed 
conditions for the piecewise functions to handle the following is a heuristic definition of 
ConnectFun. The initial step is to determine the forward transition function with boundary 
conditions u i and RG and the reverse transition function with boundary conditions RG and 
u 2 . If the length of the region is sufficiently long that each transition function can attain its 
maximum constant value, then we simply join the transition functions over the appropriate 
domains. If however, the length of the region is not long enough for each transition function 
to attain its maximum constant value, then we set the u-value equal to the minimum of the 
2 values which the distinct transition functions return. 


1.4 Finding the Trip Time 


Phew! We have finally solved for u(s). Here is an example output for a driving route from 
San Francisco to Los Angeles under the .5 g limit: 
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And this output is for a driving route from New York to Boston: 



Notice how the trip time is on the right of each image. How do we do this? Well, with our 
velocity profile v(s) in hand, we can simply calculate the trip time by an integral, letting s 
range from 0 to S, where S is the length of the route: 
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Chapter 2 

Calculating Money M(^) 


Implicit in space, separate from 7 itself is the notion of the cost m as a function of location, 
m : M 2 — > M. Then 

https : / /writelatex.s3.amazonaws.com/hncypvvzcbqq/page/f82b09d3b9be4e73b3ed9050daffc 

We throw out impossible locations by defining m piecewise: 

{ 00 if impassible 
20 if farm 
5 if highway 

For cost we divided the various components into a set that depends solely on length {/(/)} 


Pylons (2.1) 

Tube construction (2.2) 

Linear Accelerators? (2-3) 

Solar Panels (2.4) 

and a set that depends on position {/(7(s))} 

Tunnel (2.5) 

Permits and Land (2-6) 


2.1 Research 

(a) Get response from Musk as to cost of individual components and land considerations 

(b) Solve differential equation for v 
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